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Abstract  A  three-cluster  model  developed  for  ternary  fission  studies 
has  been  applied  for  the  first  time  to  study  the  three-body  structure  of  2n  and 
2p  halo  nuclei.  For  the  experimentally  known  2n,  2p  halo  nuclei,  all 
possible  ternary  fragmentation  potential  energy  surface  (PES)  is 
calculated.  The  two-body  breakup  reported  earlier,  clearly  indicated  a 
strong  minimum  in  the  PES  corresponding  to  In/ lp  and/or  2n/2p  cluster 
plus  core  configuration.  However,  the  present  calculations  of  PES  reveal  that, 
the  three-  body  breakup  does  not  result  always  with  2n  and/or  2p  as  a 
cluster.  A  In  and/or  lp  cluster  along  with  the  core  is  initially  formed,  and 
then  the  core  loses  one  nucleon  to  make  either  a  2n  plus  core  or  2p  plus  core 
structure.  The  results  are  substantiated  with  the  calculations  of  preformation 
probability  calculated  within  quantum  mechanical  fragmentation  theory. 
Keywords:  Halo  nuclei,  Cluster  Core  model,  Preformation  Probability 


1.  INTRODUCTION 


The  term  ‘halo’  is  familiar  to  all  of  us  and  is  used  in  wide  areas  with  different 
meanings.  In  nuclear  physics,  a  halo  nucleus  means,  the  nucleus  with  a  large 
spatial  extension  caused  by  protons  or  neutrons  around  the  central  core.  The  idea 
about  the  nuclear  halo  structure  first  came  into  focus  through  the 
experiments  done  by  Tanihata  et  al.  In  an  experiment  made  at  Lawrence 
Berkeley  Laboratory  to  measure  the  interaction  cross  sections  and  nuclear  radii 
of  various  nuclei,  they  identified  the  nLi  nucleus  with  a  large  radius  and 
deformation.  These  nuclei  with  larger  radius  are  termed  as  “Halo”  nuclei 
and  are  found  to  exist  at  the  extreme  edges  of  the  table  of  nuclides- v/z., 
the  neutron  and  proton  drip-line.  They  are  characterized  by  a  very  small 
separation  energy,  and  very  short  half-lives  (milliseconds).  The  separation 
energy  of  the  outer  nucleons  is  really  small  and  nuclear  halo  appears  once 
when  it  has  small  binding  energy  and  low  angular  momentum.  Two-body  halos 
occur  for  nucleons  in  s-  or  p-  states.  Three-body  halos  form  with  hyper- 
spherical  quantum  numbers  K  =  0  or  K  =  1 . 
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Sreeja,  I  After  the  experimental  reports  of  Tanihata  et  al.,  [1-3]  several  studies  of 

Balasubramaniam,  M  neutron/proton  halos  are  done  both  experimentally  and  theoretically.  Direct 
mass  measurement  of  the  two-neutron  halo  nucleus  6He  and  improved  mass 
measurement  for  the  four-neutron  halo  8He  are  done  using  the  TITAN  Penning 
trap  mass  spectrometer  at  the  ISAC  facility  which  helps  for  the  precise  test 
of  three-body  forces  at  neutron-rich  extremes  [4],  Danilin  et  al.  [5]  studied 
the  halo  structure  of  the  three-body  Borromean  continuum  theory  using  the 
method  of  hyper-harmonics  with  core  +  n  +  n.  Gupta  et  al.  [6-9]  studied  the 
halo  structure  of  neutron  and  proton  drip  line  nuclei  for  both  ln/2n  and  lp/2p 
using  the  cluster  -  core  model.  They  reported  that  the  stability  of  the  core 
nucleus  is  found  to  be  at  both  N=2Z  and  2Z  ±  2.  Three-body  halo  structure 
of  22C  (20C  +  n  +  n)  is  studied  by  Horoiuchi  et  al.  [10]  using  the  orthogonality 
constraint.  Similar  studies  are  carried  out  for  uLi  and  l4Be  using  zero-range 
interactions  [11].  Two  neutron  halo  structure  of  20C  is  studied  using  three-body 
formalism,  which  showed  the  occurrence  of  at  least  two  Efimov  states  [12]. 
Using  the  microscopic  self-consistent  relativistic  mean  field  description,  Shi- 
Sheing  Zhang  et  al.,  predicted  a  p-orbit  In  halo  structure  for  31  Ne  and  found 
that  there  is  no  evidence  of  an  s-orbit  In  halo  in  29Ne  [13].  The  halo  structure  of 
Borromean  systems  like  "Li,  14Be  and  22C  are  studied  by  Souza  et  al.  [14].  Two 
neutron  halo  structure  of  22C  is  studied  using  a  Glauber  model  [15].  Recently 
an  effective  field  theory  has  been  used  to  study  the  three-body  halo  structure 
of  6He  [16]. 

In  Ref.  [17,18],  the  concept  of  low  proton  separation  energy  and  abnormal 
occupation  of  outside  protons  is  used  for  the  studies  of  proton  halos.  Using 
the  nonlinear  RMF  theory  Ren  et  al.  [19]  studied  the  halo  structure  of  26P  (lp 
halo)  and  27S  (2p  halo).  The  first  experimental  measurement  of  2p  removal  was 
done  by  Kanungo  et  al.  [20,21]  using  I7Ne.  Grigorenko  et  al.  [22]  studied  the 
possibility  of  two-proton  halo  in  17Ne  using  three-body  models.  Geithner  et  al. 
[23]  experimentally  measured  the  charge  radii  of  18_23Ne  for  the  first  time  and 
figured  out  the  2p  halo  structure  of  17Ne  using  Penning  trap  mass  spectrometry 
and  collinear  laser  spectroscopy.  A  recent  study  [24]  of  the  2p  halo  structure 
of  17Ne  using  MCM  disagrees  with  the  previous  results. 

In  the  present  work,  encouraged  by  the  success  of  cluster-core  model,  in 
the  a  priori  prediction  of  halo  nuclei,  we  extended  the  ideas  of  CCM  into  the 
three-cluster  model  proposed  by  one  of  us  for  the  studies  of  ternary  fission. 
Within  three-cluster  model,  the  three-body  breakup  of  all  the  known  2n/2p 
halo  nuclei  are  studied.  Further,  the  charge  minimized  ternary  fragmentation 
potential  of  halo  nuclei  is  used  to  calculate  the  preformation  probability  for 
the  first  time  within  quantum  mechanical  fragmentation  theory.  In  section  2,  a 
brief  account  on  the  cluster-core  model,  and  three-cluster  model  is  discussed. 
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In  section  3,  the  obtained  results  are  presented  and  discussed.  Finally,  in 
section  4,  the  summary  and  conclusions  are  outlined. 

2.  METHOD 

The  cluster-core  model  (CCM),  of  Gupta  and  collaborators  [6-9]  is  an 
extension  of  preformed  cluster-core  model  [25-29].  Earlier  using  CCM,  the 
two  body  fragmentation  potential  energy  surfaces  (PES)  of  known  halo  nuclei 
were  studied.  Taking  clue  from  CCM,  for  the  first  time,  we  extended  the 
study  to  the  three-body  fragmentation  using  the  three-cluster  model  (TCM) 
[30-37]. 

2.1  The  cluster-core  model  (CCM) 

In  CCM,  the  binary  fragmentation  potential  between  the  two  fragments  is 
calculated  as 


v  =y2 

'  tot  ^-ti- 


<  +  vc  +  vP 


(1) 


Here  mx  are  the  mass  excesses  of  the  fragments  in  energy  units  taken 
from  Audi-Wapstra  [38]  and  Moller-Nix  [39]  tables.  For  neutron  clusters, 
the  mass  excess  is  taken  as  xAmn ,  where  x  is  the  number  of  neutrons  and 

A mn  =  8.07 1 3Mc;F.  For  proton  clusters,  it  is  taken  as  xA m  -acAJ  ,  where  x 

is  the  number  of  protons  and  A m  =  7. 2 8  80  Me  F  and  ac  =  0.7053  Me  F  . 


Vc  = 


Z{Z2e2 

R, 


(2) 


is  the  Coulomb  interaction  energy  between  the  two  fragments,  where 
Rt  =  Rt  +  R2  is  the  radius  of  the  fragments  in  the  touching  configuration, 
which  is  defined  as  Ri  =  1 .28 T2  3  -  0.76  +  0.8 A7infm  (i  -  is  taking  the  values 
of  1  and  2  for  the  two-body  fragmentation). 

According  to  Blocki  et  al.  [40]  the  proximity  potential  V  is, 

Vp  =  4t TRjb^(C).  (3) 

Here  the  mean  curvature  radius,  R  has  the  form, 
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(4) 


R  =  RlRl 


R\  +  R-, 


The  specific  nuclear  surface  tension  7  is  given  by 


7  =  0.9517 


7 

1 

N 

_ _ 1 

1-1.7826 

K.  A  )  \ 

MeVfm2 


(5) 


Here  N,  Z  and  A  are  the  neutron ,  charge  and  mass  number  of  the  parent  nucleus . 
b  is  the  diffuseness  of  the  nuclear  surface  given  by 


b  =  f  7T  /  2>/3  In  9l 

L  J  u 


(6) 


Here  t10_90  is  the  thickness  of  the  surface  in  which  the  density  profile  changes 
from  90%  to  10%.  The  value  of  b  =  0.99  fm. 


i(£ -  2.54)2  - 0.0852(5  - 2.54)’  fori  <  1 .25 1 1 


-3 .437  exp 


£ 


(7) 


0.75 


\for£  —  1-2511 


Where,  £  =  7  Here  this  function  is  zero  as  separation  distance 
b 

s{Ru  -Rl-R2)  =  0  for  touching  configuration. 

The  cluster  preformation  probability  P0  is  obtained  by  solving  the 
Schrodinger  equation  in  7  motion  at  fixed  R, 


r 


1 


2jB  dri  IB  dri 

\  nn  '  \  rm  1 


+  v{v) 


r{v)=E''r{v) 


(8) 


Solving  this  equation  numerically,  \^{ji)\  gives  the  probability  of  finding  the 
mass  fragmentation  rj  at  a  fixed  position  R  on  the  decay  path.  Normalizing 
and  scaling  to  give  the  fractional  mass  yield  at,  say,  the  mass  A2  of  the 

(  2  3 

cluster  drj  -  _  ,  the  cluster  preformation  probability,  for  the  ground  state 

l  A) 

decay,  in  this  model  is 
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(9) 


B  (r/)  is  the  mass  parameter  which  represents  the  kinetic  energy  part  in  the 
above  equations  are  the  smooth  hydrodynamical  mass  parameters  of  Kroger 
and  Scheid  [41].  This  gives  a  simple  analytical  expression  which  is  found  to 
compare  nicely  with  microscopic  cranking  model  calculations. 

For  the  B  mass  we  get, 


AmR1  v,  (1  +  0) 

~[vc(\  +  82)-\ 


with 


P  =  — 

2  R 


1 


1  +  cos6l 


(  D  ^ 

1-^ 

R 


+  - 


1 


i  / 


1  +  cosO , 


'i-A' 

R 


'2  y 


(10) 


(11) 


8  =  ^[(i ■ - cosOx)(R{  -Rc)  +  (  1  - cos02 ){R2  - Rc )]  (12) 


v  =  nR2R 

c  c 


(13) 


and  vt  =  v {  +  v2 ,  is  the  total  conserved  volume. 

The  angles  0{=02=  0  and  5  =  0,  correspond  to  two  touching  spheres. 
Rc  (V  0)  is  the  radius  of  a  cylinder  of  length  R,  having  a  homogeneous  flow  in 
it,  whose  existence  is  assumed  for  the  mass  transfer  between  the  two  spherical 
nuclei  (The  flow  in  nuclei  is  assumed  to  be  radial).  The  authors  use  a  special 
Ansatz  for  R  : 

c 


R 


c 


a.  min  yRx,  R2)f 

V 


R  ' 
0, 


(14) 


with  f(x)  =  1  for  x  <  1  and  a  =  0.4  to  0.8  for  the  best  fit  to  the  microscopic 
calculations. 
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2.2  The  three-cluster  model 

Three-body  fragmentation  potential  of  all  possible  combinations  is  calculated 
using  the  TCM  by  imposing  the  condition  Al  >  A-,  >  A3  on  the  masses  to  avoid 
repetition  of  the  combinations.  The  potential  is  calculated  using 


v  =y?  m  +  y3  y3  v 

v  tot  ^i= 1  "lx  T  ^i=  1  i=i  y  ii 


(15) 


where 


Kt=V(ii+VPii  (16) 

Z,  Z,e2 

Here,  for  the  Coulomb  interaction  energy,  VCij  =  ~ —  ,  the  radius  between 

the  fragments  R.  is  taken  as,  ^ii 


jR.^2  —  +  R~2 

R23  =R2+R3 

and 

Rn  =  R}  +  2R2  +  Rv 

And  for  the  proximity  potential,  VPj  =  Att Ripb(\,{t^)  ,  R.  has  the  form, 

-  _  Rt  R2 
12  Rt  +  R2 
—  RnR-i 

D  —  2  3 

-3  R2  +  R3 

and 

-  _  R3Rt 

31  R3  +  2  R2  +  R{ 

In  the  specific  nuclear  surface  tension  7 ,  the  terms  N,  Z  and  A  changes  for 
each  two  fragments. 
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The  universal  function  <b(£)  depends  only  on  the  surface  separation 
s  =  Rj}  -  ( Rt  +  R  j  ■  It  is  negative  in  the  overlap  region,  zero  in  the  touching 
configuration  and  is  positive  in  the  separated  configuration.  In  this 
model,  the  potential  energies  are  calculated  for  all  possible  cluster-core 
( A,  -  A, )  &  ([Aj  ,Al]-A1)  configurations  which  help  to  find  the  neutron(s)  or 
proton(s)  cluster  +  core  configuration  with  a  minimum  potential  energy. 

3.  RESULTS 

We  have  investigated  the  three-body  structure  of  In,  lp,  2n  and  2p  halo  nuclei 
within  the  three-cluster  model.  The  three-body  potential  energy  surfaces  (PES) 
calculated  for  all  the  proton-rich  and  neutron-rich  halo  nuclei  are  shown  in 
table  1 .  The  first  column  of  the  table  lists  the  parent  halo  nucleus.  The  second 
to  fourth  column  lists  the  deepest  minimum  in  the  PES  for  each  nucleus  along 
with  any  other  pronounced  minimum  of  the  cluster  1 ,  cluster  2  and  the  core 
nucleus  resulting  from  the  potential  energy  calculations.  In  column  5  and  6, 
the  Z  and  N  values  of  the  core  is  presented.  When  the  PES  was  calculated 
for  a  two-body  breakup  as  was  reported  in  [5,6],  the  PES  clearly  indicated 
a  stronger  preference  for  the  halo  cluster  plus  core  configuration.  The  halo 
configuration  corresponding  to  two  body  potential  is  listed  in  the  last  column 
of  the  table. 

If  we  consider  the  three-body  breakup  and  in  turn  the  three-body 
fragmentation  potential,  the  scenario  is  not  as  it  was  before.  For  example,  a 
2n  -  halo  like  6He  exhibits,  stronger  minima  corresponding  to  either  ln+2H+3H 
or  'n+'I  I+4I I.  This  indicates,  that,  it  forms  as  a  breakup  with  In  as  one  of  the 
clusters  and  then  In  from  either  3H  or  4H  gives  rise  to  a  2n  cluster  and  the 
remaining  2H+2H  in  the  former  case  and  1H+3H  in  the  later  case  forms  as  4He, 
finally  giving  rise  to  a  2n+4He  structure.  We  interpret  that,  2n  or  2p  halo  nuclei 
are  first  formed  as  a  three-body  structure  with  In  or  lp  cluster  and  then  would 
finally  result  in  the  amalgamation  of  remaining  two  clusters  to  form  a  core 
donating  In  and/or  lp  to  result  in  2n  and/  or  2p  cluster  plus  core. 

Figures  1  to  4  present  the  three-body  potential  energy  surface  as  defined 
in  Eq.  (15)  as  a  function  of  the  three  mass  numbers  denoted  by  A  ,  A  and  A 
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corresponding  to  some  representative  cases  of  2n  halo  ("Li  and  19  B),  In  halo 
( 1  'B  and 24 F) ,  2p  halo  (20Mg  and 27  S) ,  and  lp  halo  CB  and 26 P)  nuclei.  In  Fig .  1 ,  the 
three-body  potential  energy  surface  as  a  function  of  mass  numbers  A  , A,  and  A3 
for  the  2n-halo  nuclei  "Li  and  19 B  are  presented  in  a  four-dimensional  ternary 
contour  graph.  For  "Li,  a  stronger  minimum  in  the  potential  energy  surface 
is  seen  corresponding  to  \H  +  4 He  +  2an  .  The  three-body  breakup  results  in 
2n  as  one  of  the  cluster  and  the  remaining  cluster  and  core  will  amalgamate 
to  form  9 Li .  However,  for  19 B  two  pronounced  minima  are  seen  in  the  PES. 
The  numerically  strongest  minimum  corresponds  to  'l Li  +  4 He  +  \n  and  the 
other  minimum  corresponds  to  [Li  +  [He  +  [n  .  Since,  ]9B  is  an  established 
2n-halo,  the  first  minimum  may  be  the  preferable  path  at  the  same  time,  the 
second  possibility  in  which  either  the  least  stable  [He  could  lose  In  to  form  a 
2n  cluster  plus  core,  can  also  be  not  ruled  out.  Similar  results  are  presented  in 
figures  3  to  4  for  In,  2p  and  lp  halo  nuclei. 


Figure  1 :  Three-body  potential  energy  surface  as  a  function  of  mass  numbers 
A  ,  A,  and  A  for  the  2n-halo  nuclei  "Li  and  l9B.  Stronger  minima  in  the 
potential  energy  surface  are  labelled. 
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Figure  2:  Three-body  potential  energy  surface  as  a  function  of  mass  numbers 
A  ,  A,,  and  A,  for  the  ln-halo  nucleus  14B,  24F.  Stronger  minima  in  the  potential 
energy  surface  are  labelled. 


In  each  figure,  the  stronger  minimum  in  the  PES  is  labelled.  In  some  cases, 
more  than  one  minimum  is  present,  indicating  different  possible  rearrangements 
of  the  nucleons  finally  to  observe  them  as  either  ln/2n  or  lp/2p  halo  nuclei. 
Similar  calculations  are  done  for  all  the  nuclei  studied  and  are  listed  in  table  1 . 

To  further  substantiate  our  result,  we  present  in  Table  2,  the  preformation 
probability  values  calculated  for  some  selected  cases  of  In,  lp  and  2n  and 
2p  halo  nuclei.  The  preformation  probability  is  calculated  for  the  two-body 
fragmentation  potential  as  defined  in  Eq.  (1).  For  the  use  of  this  fragmentation 
potential  and  the  mass  parameters  defining  the  kinetic  energy  part  of  the 
Hamiltonian  as  defined  in  Eq.  (10),  the  equation  of  motion  given  in  Eq.  (8)  is 
solved.  The  square  modulus  of  the  wave  function  after  normalization  gives  the 
preformation  probability  as  denoted  in  Eq.  (9). 

It  is  clearly  seen  from  the  Table  2,  that,  ln+core  and/or  lp+core  nuclei 
are  found  to  have  a  larger  probability  of  formation  supporting  the  results  seen 
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Table  1:  The  deepest  minima  seen  in  the  three -body  potential  energy  surface  of  the 
known  In,  lp,  2n  and  2p  halo  nuclei  are  listed.  If  there  is  more  than  one  pronounced 
minimum  in  the  PES,  which  is  also  listed. 
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Cluster2 

Core 
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10 
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through  PES  for  the  two-body  breakup.  At  the  same  time,  for  2n  and/or  2p  halo 
nuclei,  the  preformation  probability  values  are  found  to  be  very  low  compared 
to  In  and/or  lp  halo  nuclei  results.  This  result  reported  for  the  first  time  for 
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Table  2:  The  preformation  probability  values  for  the  use  of  two-body  breakup  potential 
for  some  representative  cases  of  In,  lp,  2n  and  2p  halo  nuclei. 


Halo 

Parent 

Preformation  Probability 

In  -  halo 

22N 

0.932 

In  -  halo 

24  p 

0.978 

lp  -halo 

26 p 

0.932 

lp  -  halo 

27 P 

0.990 

2n-halo 

22C 

0.661 

2n-halo 

27F 

0.626 

2p-halo 

20Mg 

0.668 

2p-halo 

27S 

0.562 
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Figure  3:  Three -body  potential  energy  surface  as  a  function  of  mass  numbers  Ap 
A,  and  A,  for  the  2p  halo  nuclei  20Mg,  27S.  Stronger  minima  in  the  potential  energy 
surface  are  labelled. 
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the  halo  nuclei  indicates,  and  strengthens  our  interpretation  that,  the  2n  and  2p 
halo  structure  is  initially  formed  with  a  larger  probability  as  In  plus  core  and / 
or  lp  plus  core  and  then  the  core  nucleus,  further  looses  either  a  neutron  to 
make  2n  plus  the  core  structure  or  a  proton  to  make  a  2p  plus  core  structure. 


Figure  4:  Three-body  potential  energy  surface  as  a  function  of  mass  numbers  Aj, 
A,  and  A  ,  for  the  lp  halo  nuclei  8B  and  26P.  Stronger  minima  in  the  potential  energy 
surface  are  labelled. 


SUMMARY  AND  CONCLUSIONS 

Three-body  breakup  and  the  corresponding  three-body  fragmentation  energy  is 
calculated  for  investigating  three-body  structure  of  2n  and  2p  halo  nuclei.  Unlike 
two-body  fragmentation,  which  revealed  clearly  the  core  plus  neutron  and / 
or  proton  cluster,  three-body  PES  exhibits  more  than  one  stronger  minimum. 
Further  the  2n  /  2p  three-body  breakup  does  not  always  result  in  one  of  the 
cluster  as  2n  or  2p.  In  most  of  the  cases,  In  or  lp  cluster  is  preferred  in  the 
ternary  breakup  and  hence  the  remaining  cluster  and  the  core  nuclei  must  lose 
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